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ABSTRACT 

We apply a cross-correlation technique to infer the S > 3mJy radio luminosity 
function (RLF) from the NRAO VLA sky survey (NVSS) to z ^ 3.5. We measure 
S the over density of radio sources around spectroscopically confirmed quasars. S is 
related to the space density of radio sources at the distance of the quasars and the 
clustering strength between the two samples, hence knowledge of one constrains the 
other. Under simple assumptions we find $ oc (1 -f out to z ^ 2. Above this 

redshift the evolution slows and we constrain the evolution exponent to < 1.01 {2a). 
This behaviour is almost identical to that found by previous authors for the bright 
end of the RLF potentially indicating that we are looking at the same population. 
This suggests that the NVSS is dominated by a single population; most likely radio 
sources associated with high-excitation cold-mode accretion. Inversely, by adopting a 
previously modelled RLF we can constrain the clustering of high-redshift radio sources 
and find a clustering strength consistent with tq = 15.0 ± 2.5 Mpc up to z 3.5. This 
is inconsistent with quasars at low redshift and some measurements of the clustering 
of bright FRII sources. This behaviour is more consistent with the clustering of lower 
luminosity radio galaxies in the local universe. Our results indicate that the high- 
excitation systems dominating our sample are hosted in the most massive galaxies at 
all redshifts sampled. 

Key words: methods: statistical, galaxies: evolution, radio continuum: galaxies, 
large-scale structure of Universe,quasars: general,galaxies: active 


1 INTRODUCTION 

The standard approach to measuring the radio luminosity 
function (RLF) requires a sample with distance information 
to convert fluxes to luminosities. These distances typically 
come from a cross-match to existing optical redshift surveys. 
Millijansky radio sources have sky densities of a few 10s 
per square degree (Condon et al. 1998; Mauch et al. 2003) 
hence require wide-field (> 100 deg^) spectroscopic surveys 
to build up significant statistics in the RLF. These samples 
exist in the local Universe (e.g. 6dF, SDSS) and in combina¬ 
tion with wide-held mjy radio catalogues the local RLF has 
been shown to be a combination of two main populations: ac- 
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tive galactic nuclei (AGN) with a double power law LF (sim¬ 
ilar to that of quasars e.g. Boyle et al. 2000) at high lumi¬ 
nosities and star-forming galaxies with a Schechter function 
LF at lower luminosities (Best et al. 2005; Mauch & Sadler 
2007). 

Wide-held spectroscopic surveys are not deep enough 
to probe the overall galaxy population at higher redshift. 
Luminous red galaxies (LRGs) are bright enough to pro¬ 
duce large samples up to z ~ 0.7 (Eisenstein et al. 2001; 
Cannon et al. 2006). While this is hardly a representative 
slice of the galaxy population, local surveys show that most 
radio sources in the lO^'^ < L < 10^® W/Hz regime (that 
translates to huxes of S ~ 1 to 100 mJy at z ~ 0.7) are 
associated with massive red galaxies (Condon & Broderick 
1988; Mauch & Sadler 2007). Using LRGs, Sadler et al. 
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(2007) found evolution described well by shifting the AGN 
portion of the local RLF in the luminosity direction by 
(1 + These results were in broad agreement with 

Clewley & Jarvis (2004) who used galaxies from the Sloan 
digital sky survey to show that fainter L < 10^® WHz~^ sr~^ 
radio sources evolved more slowly than brighter ones up to 
2 ~ 0.5 

Deep pencil-beam optical surveys offer higher-redshift 
galaxy samples that, when combined with deep radio imag¬ 
ing, constrain the sub-mjy RLF. At these lower flux den¬ 
sities (< 0.1 mjy) radio surveys become dominated by 
star forming galaxies (Seymour et al. 2008; Padovani et al. 
2009) that show strong luminosity evolution oc (1 -I- ® 

(Padovani et al. 2011; McAlpine et al. 2013), with some con¬ 
tribution from radio-quiet AGN (Jarvis & Rawlings 2004; 
Simpson et al. 2006). The lower-luminosity AGN found 
in these surveys show somewhat less evolution than the 
Sadler et al. (2007) result. Padovani et al. (2011) find no 
evolution in their AGN to 2 : ~ 5, and when they remove 
possible star-formation derived emission they find negative 
evolution. They suggest this may be a result of extremely 
high redshifts objects in their sample and the RLF cutting 
off and declining for 2 > 1 — 2. Below these extreme redshifts 
Smolcic et al. (2009); McAlpine et al. (2013) find slow but 
significant evolution in their AGN: oc (l-|-z)^'^ and (l-|- 2 :)° ® 
respectively. 

Small radio samples with complete spectroscopic cov¬ 
erage constrain the bright end of the RLF at high 2 
(Dunlop & Peacock 1990; Willott et al. 2001). These stud¬ 
ies show that at bright fluxes radio sources are found up 
to high redshift ( 2 ^ 3 ), indicating the difficulty in obtain¬ 
ing complete spectroscopy on large radio samples. They also 
found strong evolution in the RLF. Willott et al. (2001) used 
a combination of tiered radio samples with a faintest limit 
of SisiMHz > 500 mJy to model the RLF. They separated 
their LF model into two populations roughly separated by 
being above or below L ~ 10^® W/Hz. The lower luminosity 
population being primarily FRI objects or FRIIs that show 
little evidence for an AGN in the optical, and the higher 
luminosity sample containing bright FRII sources often as¬ 
sociated with optical quasars. The brighter population’s LF 
increases towards higher 2 peaking at 2 ~ 2 and then falling. 
The fainter end is described by a Schechter function that 
increases until it reaches 2 ~ 1 after which it remains sta¬ 
tionary. In reality the lower luminosity population is poorly 
constrained for 2 : > 1 although further strong evolution is 
ruled out by source counts. 

Above 2 ~ 0.7 the mJy RLF is difficult to con¬ 
strain. It lies between the parameter spaces constrained by 
pencil-beam surveys that run out of radio sources at higher 
flux densities, and the targeted surveys that require large 
amounts of telescope time to push fainter. In this regime 
the RLF has been estimated from samples that have semi- 
complete spectroscopic coverage supplemented by photo¬ 
metric redshifts. Waddington et al. (2001) used a 1 mJy lim¬ 
ited sample of 72 galaxies with 65% spectroscopic complete¬ 
ness to show that the evolution of fainter ~ 10^^ W radio 
sources peaks later compared with brighter ~ 10 ^® sources. 
Rigby et al. (2011) used a tiered sample that included the 
Waddington et al. (2001) sample and pushed further down 
to 0.1 mJy at 1.4 GHz (for 2 < 1.3) using photometric red¬ 
shifts from the COSMOS field. They confirmed this differ¬ 


ential evolution with radio luminosity analogous to ‘down¬ 
sizing’ seen in star formation rates and Xray/optical AGN. 

The flux range 1 < -S' < 100 mJy is of particular in¬ 
terest since it samples the RLF in the luminosity regime 
where the bulk of the energy density from AGN is emit¬ 
ted from redshifts 0.5 < 2 < 3.5; the peak of AGN activity 
in the Universe. Hence this flux range is fundamental to 
our understanding of radio AGN and their impact on their 
surrounds. Constraining the LF in this parameter space is 
difficult and has thus far only been possible in small sam¬ 
ples with incomplete spectroscopy. In this paper we look at 
an alternative approach. We use spectroscopic quasars as a 
tracer of the large-scale structure at high redshift and cross 
correlate these with the NVSS to determine the RLF. 

Throughout this work we will assume a standard flat 
(Dm, Da) = (0.3, 0.7), h = 0.7 cosmology. The paper is or¬ 
ganised such that we give the background to our technique 
in section 2, in section 3 we introduce the data sets we will 
be using for our analysis that is described in section 4. In sec¬ 
tion 5 we show our results and discuss their meaning for the 
RLF in section 6 and high-redshift clustering in section 7. 
We summarise our results in section 8 . 


2 REVIEW OF TECHNIQUE 

The technique we follow exploits a data set that has redshift 
information to constrain a sample that does not (Phillipps 
1985; Phillipps & Shanks 1987). This approach has been 
used for many years and has been recently exploited to 
reproduce the redshift distribution of photometric samples 
(Newman 2008; Matthews & Newman 2010) and similarly 
to calibrate photometric redshifts (Schulz 2010). Here we 
briefly outline the process we will follow as described in 
Phillipps (1985) and Phillipps & Shanks (1987). 

We begin with the real-space correlation function ^(r, 2 ) 
defined such that, for a galaxy population with space density 
(j){L, 2 ), the probability of finding a galaxy in a volume 6V 
a distance r from an arbitrary galaxy is 


5P = ^{L,z)[l + aL,z,r)]SV. (1) 


In the linear haloe-haloe regime (1 r < 100 Mpc) the 
correlation function is well-described by a power law ^ = 
{ro/ry with 7 ~ 1.8 (e.g. Peebles 1980). 

The angular statistic Eeicess is defined as the excess 
number of galaxies with luminosity L to L + SL within a 
projected radius R of an arbitrary galaxy with known red¬ 
shift 2 . Assuming a power law form for the correlation func¬ 
tion, and that the evolution of ro and tj) are minimal over a 
clustering length 


-‘excess 


{L,z) 


27rG(7)ro^(L,2)R®-^0(L,2) ^^ 


( 2 ) 


where G is a constant defined by 7 (see Phillipps 1985). 

Importantly "Eexcesa is trivial to measure between a 
sample with redshifts and one without. We may then con¬ 
strain the clustering strength ro{z,L) and luminosity func¬ 
tion 0 ( 2 , L) of a population with no redshifts. 



Thp mi'll rnl'in liimirinQ'iin fu'nri'inri 


9 . 


O 



0 2 4 


z 

Figure 1. The redshift distribution of our final quasar sample. 
The fine lines show the DR7 (black) and DRIO (grey/red) samples 
separately and the heavy line shows the total. This demonstrates 
the extra redshift range from z ~ 2 to 3.5 made available by the 
inclusion of the DRIO quasars. 

3 DATA 

In this paper we aim to measure the radio luminosity func¬ 
tion and clustering strength of high-redshift radio sources. 
We do this by counting quasar-radio source pairs from a 
spectroscopic quasar sample that has redshifts and a radio 
catalogue that has none. 

We take the NVSS survey (Condon et al. 1998) as our 
parent sample of radio sources. The NVSS covers the whole 
sky north of —30°, but for our purposes we are only in¬ 
terested in extragalactic sources and so cut out all objects 
with galactic latitude |6| < 10° (as well as Dec. < —30). 
We also make a flux cut at 3mJy above which the NVSS 
is ~ 90 % complete (Condon et al. 1998) leaving 1,062,117 
radio sources in our sample, the vast majority of which have 
no distance estimate. 

The quasar sample we use is a combination of the 
Sloan Digital Sky Survey (SDSS) DR7 quasar catalogue 
(Schneider et al. 2010) and the Baryon Oscillation Spectro¬ 
scopic Survey (BOSS) DRIO quasar catalogue (Paris et al. 
2014). We combine the two samples since they cover differ¬ 
ent redshift ranges: DR7 0.1 < z < 2, DRIO 2 < z < 3.5. 
Neither sample has a single consistent selection function and 
both are somewhat unevenly distributed across the sky. To 
flatten the DR7 catalogue we follow the simple cut made by 
(Schneider et al. 2010) and only allow objects with i < 19.1: 
the magnitude limit of the main SDSS quasar survey. For 
the BOSS sample we make a similar cut of i < 20 where the 
number counts begin to turn over. Further cuts to these sam¬ 
ples are required to construct accurate random catalogues 
as discussed in the next section. 


3.1 Random catalogues 

To estimate the excess Eeaicess we need random catalogues 
matched to our radio and quasar samples to form a compar¬ 
ison. We create the random radio catalogue by generating 
random sky positions with Dec. > —30° and |b| > 10°. We 
also assign each source a flux drawn from the NVSS at ran¬ 



Figure 2. The distribution of our final quasar samples in the sky. 
Red and Blue are the SDSS DR7 and BOSS samples respectively 
while the black line shows the dec.> —30° and galactic |6| > 10° 
cuts that define the NVSS area we consider. 

dom. In case of any variation in the flux distribution of the 
NVSS due to the changing beam with declination we discre- 
tise the NVSS catalogue into one degree declination strips 
and only draw a flux value for our random source from ob¬ 
jects within the same declination strip. 

To create the random quasar catalogue we use MAN¬ 
GLE and the ‘SDSS DR72’ radial selection function from the 
mangle website (Hamilton & Tegmark 2004; Blanton et al. 
2005; Swanson et al. 2008). Note the DR72 mask was de¬ 
veloped to reproduce the sky coverage of the main DR7 
spectroscopic galaxy survey, and does not include additional 
fields that were in the DR7 quasar catalogue. Therefore we 
apply this mask to cut the area of our real DR7 quasar cata¬ 
logue as well. We then produce a random catalogue from the 
DR72 mask with ten times the number of random objects 
as quasars. 

To make the random BOSS DRIO sample we again use 
MANGLE and the same DR72 mask. Note this mask does 
not include approximately a quarter of the BOSS survey 
that was only covered photometrically after DR7. However, 
the DR72 mask reproduces the small-scale coverage of the 
survey and so we accept this loss of objects. 

To cut the DR72 area to just that observed by BOSS we 
take the field centers of the spectroscopic observations from 
the SDSS website and only include objects within 1.49° of a 
field center. Again, we produce this random catalogue with 
ten times the number of objects as BOSS quasars. 

Our final quasar catalogue has 80,494 objects, 63,682 
from DR7 and 16,812 from DRIO. Figure 3 shows the red¬ 
shift distributions of the final samples split by their survey. 
Clearly the inclusion of the BOSS DRIO quasars extends the 
redshift coverage of our sample from z ~ 2.2 to 3.5. The dis¬ 
tribution of these quasars on the sky along with the NVSS 
boundaries are shown in Figure 3 


4 ANALYSIS 

The excess number of radio sources around quasars at a 
given redshift and radio luminosity (calculated assuming 
the redshift of the quasar), T,qii{z, L), constrains the cross¬ 
clustering strength roQR{z, L) and the radio luminosity func¬ 
tion (f>{z,L) (Equation 2). By assuming prior knowledge of 
either roQR or (j) we can then constrain the other. In this 
section we describe models we will assume for toqr and our 
method for estimating E. 
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Figure 3. The clustering strengths (ro) for the populations 
we are considering. Black points give the quasar autocorrelation 
strength from Ross et al. (2009) and the solid line is our simple fit 
to these. The dashed black line gives the constant radio galaxy au¬ 
tocorrelation strength from Fine et al. (2011). Grey (red online) 
lines show the correlations strengths assumed in the SKADS sim¬ 
ulation: solid radio-quiet quasars, dashed FRIs and dotted FRIIs. 

4.1 The clustering strength of quasars and radio 
sources 

In reality the cross-clustering strength roQH that appears in 
Equation 2 is rarely measured. More commonly the autocor¬ 
relation strengths of quasars, r^qq, or radio sources, roriK, 
are studied. We will relate these quantities by assuming lin¬ 
ear bias such that r^qj^ ~ r^qq x tqrr (e.g. Wake et al. 
2008) and model roqq and torr as a function of redshift 
and luminosity based on recent analyses. 

Studies of quasar clustering have repeatedly shown 
that the quasar correlation function is roughly independent 
of quasar luminosity (da Angela et al. 2008; Shanks et al. 
2011 ). roqq increases slowly with redshift and, since the 
mass clustering is falling as redshift increases, the quasar 
bias rises quickly with redshift. Converting bias into mass 
via Press-Schechter theory implies that the average dark 
haloe mass of quasars is roughly constant with Mdh ~ 
10^^ Mq at all redshifts (Groom et al. 2005; Myers et al. 
2006; Ross et al. 2009). 

Assuming there is no variation in quasar clustering 
strength with luminosity we only need the variation with 
redshift. Figure 4.1 shows the autocorrelation clustering 
strength of quasars as a function of redshift from Ross et al. 
(2009). We perform a minimisation for evolution of the 
clustering strength assuming the quadratic form 

ro = a + bz^, (3) 

where ro has units of Mpc assuming h = 0.7. We find a = 
6.8 ± 0.31 and b = 0.63 ± 0.26. We will use this empirical fit 
to estimate roqq(z). 

The clustering of mjy radio sources has been exten¬ 
sively studies at redshifts below ~ 0.8 with samples cross- 
matched to optical spectroscopic or photometric galaxies 
(Peacock & Nicholson 1991; Brand et al. 2005; Wake et al. 
2008; Donoso et al. 2010; Fine et al. 2011; Lindsay et al. 
2014). At the radio luminosities sampled in those surveys 


{L > 10^^ W/Hz) the radio population is dominated by 
AGN typically hosted by LRGs. Fine et al. (2011) showed 
little evolution in the clustering strength of these objects 
with a constant roRR ~ 11.5Mpc (Fine et al. 2011). 
This broadly matches the lack of clustering evolution seen in 
optically selected LRGs (Bell et al. 2004; Wake et al. 2006; 
Brown et al. 2007). At redshifts greater than 0.8 there are 
few indications of the clustering of radio sources due to 
the lack of wide-held optical galaxy samples in this redshift 
regime. Lindsay et al. (2014) used photometric galaxy sam¬ 
ple over a small field with deep radio observations to measure 
the cross correlation with IR galaxies. They found no varia¬ 
tion in clustering strength with radio power and while their 
correlation increased with redshift this is primarily driven by 
evolution in their IR galaxy sample rather than the radio. 

The dependence of clustering strength on radio lumi¬ 
nosity is not well described. In their angular correlation 
analysis of the NVSS Overzier et al. (2003) derived a cor¬ 
relation scale length of roRR ^ 6 /i“^Mpc for lower lumi¬ 
nosity sources < 10^® W/Hz while the brighter, potentially 
FRII, sources had a scale length roRR ~ 14/i“^Mpc. On 
the other hand clustering analyses of radio surveys matched 
to optical galaxies have found no luminosity dependence 
in the large-scale haloe-haloe regime (Donoso et al. 2010; 
Fine et al. 2011; Lindsay et al. 2014). Note that since the 
Overzier et al. (2003) had no redshift information a con¬ 
siderable series of assumptions about the radio population 
were required to derive their result, on the other hand 
both Donoso et al. (2010) and Fine et al. (2011) struggled 
for sources in their samples with L > 10^® W/Hz while 
Lindsay et al. (2014) had none. 

Given the few constraints on roRR we initially make 
the simplest empirical assumption. That is the clustering 
strength of our radio sample is constant with redshift and 
luminosity with roRR ~ 11.5 h~^ Mpc. Assuming linear bias 
and roqq{z) from Ross et al. (2009) this makes up our em¬ 
pirical (EMP) model for roqR{L, z). 

As an alternative and check we also consider the values 
assumed in Wilman et al. (2008) when they were attempt¬ 
ing to model the radio sky. They followed Overzier et al. 
(2003) and assumed considerably stronger clustering for the 
brightest radio sources. For z < 1.5 they assumed constant 
dark haloe masses of 10^® and 10 ®"*Mq for FRI and FRII 
sources respectively. At high redshift the clustering strength 
of their FRII sources would become unphysically large and 
so for 2 > 1.5 they held the bias of their FRI and FRII 
sources constant. We make the simplistic assumption that 
all radio sources with L < 10^® W/Hz are FRI sources, the 
rest being FRIIs. For radio-quiet AGN, essentially quasars, 
they assumed a constant haloe mass of 3 x lO®^fi“®M 0 with 
a similar redshift cut at 2 = 3 above which the bias was 
held constant. We will refer to this alternative model for 
roQR as the W08 model (see Figure 4.1 for a comparison of 
the differing models). 

4.2 Removing radio-loud quasars 

The statistic Tlexcess defined in section 2 is the excess num¬ 
ber of radio sources around quasars. In the derivation of 
Equation 2 it is assumed that this excess comes only from 
the clustering of matter. Radio-loud quasars in our sample 
increase the measured E and bias our results. To illustrate 
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Figure 4. The quasar—radio source angular correlation function 
of our sample. The vertical lines are at 2 and 20 Mpc projected 
distance. The redshift limits for each bin are given in the bottom 
left of the panels. 


this effect Figure 4.1 shows the angular cross correlation 
function wqr{6, z) for our sample of quasars (split into red- 
shift bins) and the NVSS catalogue. The vertical lines in 
the figure show 2 and 20 Mpc projected onto the angular 
scale at the redshift of the bins. The upturn below ~ 2 Mpc 
is caused by a combination of radio-loud quasars and non¬ 
linear single haloe clustering (e.g. Blake & Wall 2002). We 
remove the radio-loud contribution by only counting pairs 
in the annulus between R = 2 and 20 Mpc. In this region the 
angular correlation function is well fit by a single power law 
indicating that the spatial correlation function ^(r) is also 
approximately a power law. We choose 2 Mpc as the lower 
limit both from inspection of Figure 4.1 and since this cor¬ 
responds to roughly the largest known giant radio galaxies 
(Saripalli et al. 2005). 


4.3 Calculating E 

To estimate E we count all radio sources with a projected 
distance between 2 and 20 Mpc from a quasar in our samples, 
Ndd- In addition to these data-data pairs we also substitute 
our random catalogues and count Ndr, Nrd and Nrr. This 
is done in redshift and luminosity bins where the luminosity 
of the radio sources are calculated assuming the redshift of 
the quasar. Redshift bins are equally spaced over the interval 
sampled by our quasars 0.1 < 2 < 3.6. Luminosity bins are 
logarithmically spaced over three orders of magnitude, the 
lower limit of which is the lowest luminosity observable in 
that redshift bin. The flux limit of the NVSS catalogue gives 
a Malmquist bias. Hence in our summations each pair is 
weighted by Vun j^max • 

Following Hamilton (1993) we estimate S with 

E = -z^{Ndd — NrrNdr/Nrd), (4) 

where Nq is the total number of quasars in the redshift bin. 
From Equation 2 we relate E to the luminosity function and 
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Figure 5. The average excess of radio sources around quasars, 
E, binned by redshift and luminosity. We show the plot with 
both a linear and logarithmic scale since several of the points are 
scattered below E = 0 by noise. 

clustering strength with 

^ 2nG{-f)r2{L,z){R^r^T. - Rf;ML,z) ,, 

^excessZ) — -- OIj. 

3-7 

(5) 

Throughout this paper we will assume 7 = 1.8 hence G — 
3.678, and use Rmax,Rmin = 20, 2 Mpc. Equation 5 becomes 

Eeaicess (L, z) = 657roQH(I/, z)cj>{L, z)5L. (6) 

We assume this simple relationship between E, ro and 0 
throughout the rest of this work. 

We use jack knife resampling to estimate our errors 
by splitting our sample into 20 even sized (by number of 
quasars) sub fields by right ascension. We then calculate E 
in each sub field and estimate the ‘field-to-field’ errors from 
the rms of these values for E (e.g. Sawangwit et al. 2011; 
Fine et al. 2011). 


5 EXCESS PAIR COUNTS 

Figure 5 shows the values of T^^xcess we calculate from our 
sample for six redshift and five luminosity bins. The way E is 
calculated means that it can be scattered to negative values 
due to noise. Hence we show both a linear and logarithmic 
scale to illustrate how the measured values and their errors 
behave. The points with E < 0 and their errors still contain 
information about our sample and need to be included in any 
analysis to avoid introducing bias. Furthermore the error 
bars are symmetric and approximately Gaussian in linear 
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space. Hence, while plots may be in log space, any fitting to 
the data is performed in linear space. 

It is apparent from Figure 5 that at hxed radio lumi¬ 
nosity E(L, z) increases slightly with redshift. This indicates 
that one of rocjij or (p is increasing with redshift. 


6 THE RADIO LUMINOSITY FUNCTION 

Figure 5 shows the radio luminosity function we calculate 
for six redshift bins. We show the LF calculated assuming 
the EMP (solid points) and W08 (open points) clustering 
models. There is very little overall difference in the LFs be¬ 
tween the clustering models. The solid lines in Figure 5 show 
the Willott et al. (2001) RTF model at the mid-point red¬ 
shift of the bin and it is clear that in general our results are 
consistent with their model. 

To describe our data further we initially fit an evolving 
power law (p = A{1 + z)“(L/10^®)^ and find a = 1.00± 0.35 
for the EMP model and 1.02 ± 1.05 for W08. However, we 
find the large redshift and flux range we sample mean this is 
not an accurate model for our data. To better illustrate the 
redshift evolution in our data we fix the luminosity exponent 
to that from our fit (/3 = —0.99EMP; -0.85 W08) and just 
£t for the amplitude of the power law in each redshift bin. 
Figure 6.1 shows the amplitude of the fitted power law at 
10^® W/Hz as a function of redshift. 


6.1 Redshift cutoff 


Figure 6.1 indicates that the increase in space density slows 
and may turnover at higher redshifts. To include this be¬ 
haviour we introduce a redshift limit and separate parameter 
for high-redshift evolution 


HL,z) = 


(l-bz)“'(L/102®)^ 

(1 + 2lim)“'"“'*(l + 


z s; Zlim 


Z > Zlim. 

(7) 


Where ah and o; are the evolution parameters above and 
below ziim- Since there can be relatively rapid evolution we 
bin our data into 25 redshift and 10 luminosity bins. We fit 
our model with simple Markov chain Monte Carlo (MCMC) 
routine iterated 500,000 times to find the preferred values. 
We £x /I = — 1 in the fitting since there can be a degen¬ 
eracy between the evolution parameter and /3 due to our 
LFs being defined in different parts of luminosity space at 
different redshifts. Figure 6.1 shows the probability distri¬ 
butions from our fitting using the EMP model, along with 
the best fit values. Clearly our data support strong a;, ~ 4 
evolution up to zum ~ 2. At high redshift we only have 
an upper limit on the evolution parameter but can show 
at least that the increase in space density stops, or turns 
over. Interestingly this redshift cut-off is almost identical to 
that found by Willott et al. (2001) for their high-luminosity 
objects Zcut-off = 1.91 ±0.16. 

To compare further with Willott et al. (2001) we fit¬ 
ted out data using a parametrisation based on their models. 
Their model ‘C’ is split into to population roughly separated 
at L = 26 W/Hz. We fit for (1 ± z)“ evolution with a red¬ 
shift cut as in Equation 7 for each population. In this fitting 
we convert from their cosmology to our own to make the 
fitted values comparable. In the high luminosity regime we 
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Figure 7. The evolution of the amplitude of the RLF from 
2 : ~ 3.5. The solid points are using the EMP model, the open 
points W08. 



A = - 7.24 ± 0.26 
aj = 3.71 ± 0.74 

ah < 1.01 (2a) 


^lim 


= 1.95 ± 0.22 
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Figure 8. The results of our MCMC fitting of our RLF data. 
Blue and red lines show the 1 and 2a marginalised constrains 
respectively. We constrain all but the high redshift evolution pa¬ 
rameter for which we can only obtain an upper limit.. 


have very few quasar-radio source pairs and consequently 
this part of parameter space is poorly constrained. On the 
other hand at low luminosities we find (i_|_2j2.35±o.89 
tion to redshift 1.94T0.43 above which we can only estimate 
an upper limit for the evolution parameter an < 1.3 (2(t). 
This contrasts with their findings of ai = 3.5 up to a red¬ 
shift cutoff at 0.7. Willott et al. (2001) have a considerably 
brighter sample than we use here and their redshift cut is 
imposed by their flux limit. It may be that since we are able 
to better define z^ut this explains the smaller discrepancy 
between our evolution parameters. 
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Figure 6. The radio luminosity function. Solid points assume the EMP clustering model while open points are W08. The solid line 
shows the Willott et al. (2001) RLF model and the dashed line shows our evolving power law model fit to the data at the midpoints of 
the bin (top-right in each panel). The redshift range for each bin is given at the top-right of each panel. The poor fit in the first redshift 
bin is due to the data being dominated by radio sources at the high-z limit of the bin. 


6.2 Discussion 

Our results are consistent with a model that evolves strongly, 
$ oc (1 -I- to t; ~ 2 above which the LF either stays 

constant or falls. Interestingly, this is approximately the 
same redshift evolution that Willott et al. (2001) found for 
their high-luminosity population of sources. The indication 
is that, rather than having two separate populations with a 
transition at L ~ 10^® W, the radio LF may be dominated 
by a single population in the luminosity-redshift regime we 
are sampling. 

Recent studies of the RLF have focused on the ac¬ 
cretion mechanisms that launch the radio jet and the 
role that the AGN may play in heating the intergalactic 
medium. Terminologies differ but we will refer to high- 
excitation radio galaxies (HERGs) associated with high- 
accretion rate optical AGN and low-excitation radio galax¬ 
ies (LERGs) associated with substantially lower accretion 
rates via advection-dominated accretion flows in massive el¬ 
liptical galaxies. At low redshifts LERGs dominate the LF 
below L ~ 10^®W (Hardcastle et al. 2007; Best & Heckman 
2012). At these fainter luminosities and lower redshifts the 
LF has been shown to only evolve slowly (Clewley & Jarvis 
2004; Sadler et al. 2007; Smolcic et al. 2009; McAlpine et al. 


2013) However, there is evidence that the HERG popula¬ 
tion evolves considerably more strongly than the LERGs 
(Willott et al. 2001; Best et al. 2014), potentially becoming 
the dominant population in the luminosity range we sample 
around z ^ 1. 

The 3 mJy flux limit we impose allows us to sample lu¬ 
minosities of 10^® W up to z = 2.5. However, the strong 
evolution of the LF coupled with the increased comoving 
volume at high redshift means that we are dominated by 
sources at z ~ 2. Assuming the simplistic power law LF 
from our MGMC fit, less than 5 % of our sample is at z < 1. 
The indication is that our signal is dominated by HERGs, 
and hence it may be unsurprising that we find almost iden¬ 
tical evolution parameters to the factor of ~ 50 brighter 
Willott et al. (2001) sample. 

Waddington et al. (2001) and Rigby et al. (2011) found 
that the fainter end of the RLF peaked in density at lower 
redshift. Due to the nature of our analysis we are always 
dominated by the radio sources close to our flux limit. Hence 
we cannot split the sample into luminosity bins to compare 
across a range of redshifts. We find a redshift cutoff ai z = 
1.95 ± 0.22. At this redshift our 3mJy flux limit translates 
to log(L/W) = 25.77, and so we can consider the turnover 
seen in our data to be due to radio sources at or somewhat 
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Figure 9. The points with error bars show tqqh as a function 
of redshift for our sample. The dashed black line gives the im¬ 
plied value of ronn assuming the empirical values of tqqq from 
(Ross et al. 2009). Grey (red online) dashed and dotted lines show 
our EMP model values for tqqq and rgnn respectively. 

brighter than this. At these luminosities Rigby et al. (2011) 
found a redshift cut closer to 2 = 1 although our results are 
consistent within a few sigma. 


7 THE CLUSTERING OF RADIO SOURCES 

Reversing what we have done above we can integrate the 
Willott et al. (2001) luminosity function above our flux limit 
to give 4> in Equation 6 and hence constrain the clustering 
strength between the radio sources and quasars in our sam¬ 
ple. Figure 7 shows the measured cross correlation strength 
in six redshift bins. While there may be some hint of an in¬ 
crease in clustering strength with z, toqr is consistent with 
a constant value of 10.4±2.5 Mpc over the full redshift range 
sampled (if poorly constrained at the highest redshifts). The 
dashed line in Figure 7 shows the value of vorr assuming 
the empirical ht in Figure 4.1 and Tqqr — roggroRR. Again 
the estimates for vqrr are consistent with a constant value 
of 15.4. 

At lower redshift {z < 1.5) we find considerably stronger 
clustering in our sample compared to quasars, more in line 
with the results for radio galaxies from Fine et al. (2011) and 
Lindsay et al. (2014). At higher redshifts both our errors and 
the clustering strength of quasars increase and we cannot 
form a distinction. 

Our results for tqrr are not consistent with the strong 
20 — 25 Mpc values assumed in Wilman et al. (2008) for 
FRII sources. This is despite our being dominated by bright 
L > 10^® W sources at high redshift and our sample po¬ 
tentially being dominated by HERGs/FRIIs at all redshift 
as discussed in our RLF analysis. None the less, we find 
are strong enough clustering to indicate these radio sources 
are in some of the most massive haloes at all redshifts we 
sample (~ 10^^ Mq at 2 ~ 0 to ~ 10^^ ® Mq at 2 3). A 

possible explanation for this would be a later (2 < 1.5) break 
imposed in the Wilman et al. (2008) bias/clustering model, 
combined with our being dominated by fainter FRI/LERG 
sources at low redshift. Alternatively the clustering strength 


could depend strongly on luminosity and redshift to contrive 
to give our results, although this has not been noted before. 


8 SUMMARY 

We measure the overdensity E( 2 , L) of radio sources around 
spectroscpic quasars and relate this to evolution in the radio 
source population from 2 3.5 to today. Our key results can 

be summarised: 

E( 2 , L) is measured in redshift/luminosity bins 
and we find significant evolution with redshift. This 
can only be explained by either the RLF of clustering 
strength increasing to 2 2. 

Under some simple models for ro{z,L) we find 
strong evolution, ()) oc (1 -|- 2 )®'^^®'^, up to 2 = 1.9 ± 0.2 
above which the evolution declines, although we can 
only constrain an upper limit. These evolution parame¬ 
ters are consistent with those found by Willott et al. (2001) 
for the brighter radio source population. The indication may 
be that the same population of HERGs dominates the NVSS 
at all flux densities above 2 ~ 1. 

Assuming the Willott et al. (2001) LF model we 
find the clustering strength of radio sources to be 
consistent with a value of tqrr — 15.0 ± 2.5 This is 
inconsistent with quasars at low redshift and the W08 model 
for FRII clustering at intermediate (1 < 2 < 2). A possible 
explanation would be the population being dominated by 
LERGs at low redshift and clustering more like quasars at 
higher redshift. Regardless, our results show that these radio 
sources are found in the most massive dark matter haloes at 
all redshift we sample. 

In this work we have demonstrated a technique that 
exploits a well defined sample with distance information to 
constrain the luminosity function and clustering of a sample 
without. The next generation of radio surveys will push still 
deeper beyond the flux limits of the NVSS used here. Despite 
new wide-held redshift surveys (e.g. EUGLID) the vast ma¬ 
jority of the sources detected in these surveys will not have 
reliable distances. The method presented in this work offers 
an alternative approach to studying these populations with 
observational strategies that are already possible and, for 
the most part, have already been carried out. 
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